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We consider resonant scatterers with large scattering cross-sections in graphene that are produced 
by a gated disk or a vacancy, and show that a gated ring can be engineered to produce an efficient 
electron cloak. We also demonstrate that this same scheme can be applied to tune the direction of 
electron flow. Our analysis is based on a partial-wave expansion of the electronic wave-functions in 
the continuum approximation, described by the Dirac equation. Using a symmetrized version of the 
massless Dirac equation, we derive a general condition for the cloaking of a scatterer by a potential 
with radial symmetry. We also perform tight-binding calculations to show that our findings are 
robust against the presence of disorder in the gate potential. 


I. INTRODUCTION 

Analogies between light and electrons are known to be 
at the origin of a plethora of interesting physical phe¬ 
nomena in condensed matter and in optics. Interference 
effects such as Anderson localization [T] and Fano reso¬ 
nances [2] are just two examples of effects that have been 
observed in both areas of Physics. In the last decade the 
advent of metamaterials and the isolation graphene have 
been fostering the discovery of novel fundamental phe¬ 
nomena and applications based on the analogy between 
light and electrons. 

On one hand, the low-energy electronic excitations of 
graphene [3], are described by Dirac fermions with linear 
dispersion relation similar to photons [4j. This “optics¬ 
like” dispersion relation has been exploited in the devel¬ 
opment of electronic analogues of photonic devices such 
as beam-splitters, wave-guides, and Fabry-Perot inter¬ 
ferometers in ballistic graphene Bm- Graphene-based 
electronic analogues of optical applications grounded on 
metamaterials, such as perfect lenses, have been also pro¬ 
posed [ilIH]. 

On the other hand progresses in the field of meta¬ 
materials [H] have allowed the development of elec¬ 
tromagnetic cloaks to achieve invisibility, and stimu¬ 
lated the quest for their electronic analogues. Among 
the existing approaches to achieve invisibility one can 
highlight the coordinate-transformation method HDHH] 
and the scattering cancellation technique p^U20] . The 
coordinate-transformation method requires metamateri¬ 
als with anisotropic and inhomogeneous profiles, which 
are able to bend the incoming electromagnetic radiation 
around a given region of space, rendering it invisible to 
an external observer. This method has been first real¬ 
ized experimentally for microwaves [12], and later for in¬ 
frared and visible frequencies mm- Extensions of this 
method to matter waves do exist [23]-[26|, proving that 
it is possible to reroute the probability flow around a 
scatterer. However, as in the electromagnetic case these 
electronic analogues require cloaking layers with extreme 
material parameters, such as anisotropic potential pro¬ 


files and/or infinite values for the effective mass, impos¬ 
ing limitations to practical implementations in condensed 
matter systems. Alternatively, the scattering cancella¬ 
tion technique for electromagnetic waves relaxes many of 
these constraints in terms of metamaterial requirements. 
It uses an homogeneous and isotropic plasmonic layer 
in which the induced polarization vector is out-of-phase 
with respect to the electric field so that the in-phase con¬ 
tribution given by the cloaked object may be partially or 
totally canceled [m m na mi. The scattering cancela¬ 
tion technique was first experimentally implemented for 
microwaves [laiiH], and later extended to acoustic [28], 
and matter waves [29]-[3T] . One of the proposals for mat¬ 
ter waves cloak, based on the expansion of partial waves 
and later extended to graphene [32], demonstrates that 
core-shell nanoparticles embedded in a host semiconduc¬ 
tor with a size comparable to electronic wavelengths can 
be “invisible” to incoming electrons [29]. To achieve this 
effect the contribution of the first two partial waves must 
vanish simultaneously so that the total scattering cross 
section can be smaller than 0.01% of the physical cross 
section. As we shall see, a similar effect can be exploited 
in graphene for the cloaking of electronic scatterers. 

Charge carriers in graphene have high mobility, and 
their peculiar properties cause small transport cross- 
sections as a result of back-scattering suppression. How¬ 
ever, different types of defects, such as vacancies and co¬ 
valently bonded adatoms, can form resonant states at, or 
close to, the Dirac point [33H35] . Resonant states increase 
the scattering cross-section by many orders of magnitude 
and hence may limit the carrier’s mobilities. They can 
also can be used to modify the electronic properties in 
the vicinity of the Dirac point. 

In the present article we combine the idea of “electronic 
cloaking” with graphene’s peculiar electronic properties 
to propose an alternative and efficient scheme for achiev¬ 
ing electronic invisibility of resonant scatterers. We show 
that sharp transitions between invisible and resonant 
states occur for some values of the potential, which can 
be viewed as an electronic analogue of optical comblike 
states in coated scatterers [36]. We consider two dif- 
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ferent resonant scatterers: an engineered resonant state 
produced by a top gate and a vacancy. We demonstrate 
that the “electron cloak” scheme proposed in Refs. [^[32] 
is strongly facilitated by the fact that the first two phase 
shifts, which give the dominant contribution in the treat¬ 
ment of electron scattering by a core-shell impurity, are 
identical for graphene. As a result, one just need to tune 
one of the scattering potentials to vanish the scattering 
cross-section, in contrast to the case where the impu¬ 
rity is embedded in a host semiconductor [29] . When the 
sub-lattice symmetry is broken in graphene, e.g. as in the 
case a vacancy, the first two phase shifts are no longer 
equal so that energy selective cloaking is more difficult 
to achieve. However, it is highly efficient at the Dirac 
point, where the resonant state is located. We also de¬ 
rive a general condition for the cloaking of a scatterer by 
a potential with radial symmetry using a symmetrized 
version of the massless Dirac equation. We extend the 
analysis to tight-binding calculations and show that our 
results are robust against the presence of disorder. Our 
findings reveal that electron cloaks are not only more eas¬ 
ily implemented in graphene, but also that they could be 
applied to selectively tune the direction of electron flow 
in this material. We suggest that new devices may be 
engineered by exploiting the cloaking effect. 

The article is organized as follows: in section |T^ we 
introduce the partial-waves formalism for graphene, and 
find the conditions for resonant scattering produced by a 
gate potential and a vacancy. In section |m] we analyze 
the effect of surrounding the two types of resonant scat¬ 
terers with a potential ring, leading to cloaking. We also 
derive a general condition for cloaking by a potential of 
radial symmetry; the example of cloaking of a void by a 
(i—function potential is provided. In section [TV| we per¬ 
form tight-binding calculations to demonstrate that our 
findings are robust against disorder in the ring, even in 
the presence of non-perfect potentials. 


II. SINGLE POTENTIAL AND THE 
RESONANT SCATTERING REGIME 

We consider a monolayer of graphene with a gated disk 
that works as a scattering center, and we choose its po¬ 
tential to maximize its total cross section. This can be 
achieved in the resonant scattering regime. Our starting 
point is the continuum-limit Hamiltonian of graphene 

no = hVEir^a^Px + (^vPy), ( 1 ) 

where p = (px^Py) is the momentum operator around 
one of the two inequivalent Dirac points K and K\ 
vf ^ 10^ m/s is the Fermi velocity, a and r denote Pauli 
matrices, with = ±1 (r^ = ±1) describing states on 
A-B sublattice (at K-K'). For such large scatterers inter¬ 
valley scattering is negligible and, in the long wavelength 
limit, for potentials with radial symmetry, the scatterer is 


described by Hy = V&{R—r) where H is a constant, R is 
the radius of the scatterer and 0(.) is the Heaviside func¬ 
tion. Here we use the partial-wave method to calculate 
the cross section of scattered waves by a radially sym¬ 
metric potential. In what follows, we derive the partial- 
wave scattering amplitudes and from their knowledge, all 
gauge-invariant quantities can be determined unequivo¬ 
cally. In cylindrical coordinates, the components of the 
graphene spinor T(r) = ('0^(r),'0^(r))^ are decomposed 
in terms of radial harmonics 

oo 

i>A{r)= E (2) 

m= —OO 

oo 

^B(r)= E (3) 

m= —OO 


where 0 = arg(/ca, + iky), k is the wavevector, m is the 
angular momentum quantum number and A{B) repre¬ 
sents the sublattice. After separating the variables of the 
graphene plus impurity Hamiltonian R/ -1- Hy, we obtain 
two coupled first order equations for the radial functions 
g^{r) and g^{r), where are Pauli matrices for sub¬ 
lattice and valley respectively. Tz is conserved and thus 
we can focus on states at K valley; scattering ampli¬ 
tudes for states at K' are quantitatively the same. In 
a partial-waves expansion, the asymptotic form for the 
spinor wave-function is given by [31 EH 

V’A,k(r) = (^ + A ) e*''’- 

(4) 

where A = ±1 denotes the carrier polarity and f{0) is 
scattering amplitude. Using a partial wave analysis, we 
can relate f{0) with the phase-shifts 6m- 

nr 

= ^ (5) 

m = — rvn 


The differential cross section is given by ^ = \f{0)\‘^ 
and the total cross section a{kR) and transport cross- 
section (jT{kR) are given by 


/ 7 A 

m=—oo 

-^[1 - cos{9)]d9 = - E sin^ ((^m - (^m+i). 


m= — oo 


(7) 


The d. c. conductivity of graphene layer at zero temper¬ 
ature can be obtained directly from arikR) |34| . 

To identify the phase-shifts, we write the spinors for 
the region inside and outside the potential as a super¬ 
position of angular harmonics. For r > R, we have 
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E = Xhvpkout and the partial-wave m is a sum of an 
incoming and a scattered wave according to 


rr:\r,e) = A^ 

+ Bm 


iAoutJm+i(fcoutr)e*(-+i)<' ) 
aoutlA„+i(fcoutr)e*(™+i)<' ) ’ 


(8) 


whereas for r < i? we have 


C(r,0)=C„ 




(9) 


(a) 
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where = \E — V\/hvF> The boundary condition 
6>) = 6>) leads to two equations. The so¬ 

lution of these equations determine the ratio 
It is straightforward to show that the phase-shift 5m 
for partial-wave m relates to the Bmj according to 
Bml^m = — tan((5^) (for more details about partial- 
waves in graphene, see references [211371 EH]). 

The Dirac equation with a radial scalar potential is 
symmetric under exchanging and which leads 

to the absence of back-scattering. The latter also corre¬ 
sponds to the relation S-m = ^m+i between phase-shifts. 
For small energies kR 1, channels m = —1,0 give the 
main contributions to f{ 0 ). 

For a fixed energy of the incoming electron, we can tune 
the potential V to produce a resonant scatterer that will 
effectively trap the electrons inside the disk, which occurs 
for phase-shifts of ±7r/2. Figure (a) shows the phase 
shifts of three partial waves; there are two resonances at 
different values of electrostatic potential, associated to 
the contributions of (io = and 6 i = 5 - 2 - 

The transport cross-section as a function of the po¬ 
tential (Figure (b)) presents pronounced peaks that 
occur exactly at the maximum of the phase- shifts, char¬ 
acterizing the resonant scattering regime. Figure [^b) 
reveals that the transport cross-section exhibits asym¬ 
metric, Fano-like resonances associated to the modes for 
m = —2,1. 
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Figure 1: (a) phase shifts (5o (black), (red) and 62 
(blue) and (b) transport cross-section (Jt ( 6 first partial 
waves are considered) as a function of the electrostatic 
potential 17 for = 0.02 eV and R = 10 nm. 


this case, if we apply the boundary conditions to Eq. 
we obtain 


6 m = — arctan 


\^m(^out-^) / 


( 10 ) 


where E = Aout^'^F^out- The zig-zag edges break sub¬ 
lattice symmetry and the relation between the Sm and 
is not valid anymore. As it can be seen from Fig.[^ 
the main contribution to the scattering cross-section at 
low energies is given by (io [38], that gives rise to resonant 
scattering for E —> 0. 


A. Circular void : another source of resonant 
scattering 


III. CLOAKING RESONANT SCATTERERS 


There are various ways to induce resonant scattering 
in graphene. One possibility is the creation of vacancies 
that give rise to zero-mode states [33l [39l [40] . Vacancies 
can be engineered by irradiating high-energy ions mm- 
This method allows the controlled incorporation of reso¬ 
nant scatterers, and can be used to taylor the electrical, 
mechanical, and magnetic properties of graphene [311142]. 
A vacancy in graphene can be modeled by a void with 
zig-zag edges [38l [43] . The boundary conditions for the 
wave-function at a void of radius R with zig-zag edges is 
given by '0 a(^ = i?) = 0 for the A spinor component. In 


After preparing the system in the resonant scattering 
regime, which maximize the transport cross-section, we 
now proceed to discuss the cloaking scheme for the scalar 
potential. We keep the original potential V (now called 
Vi) in the disk of radius R (now called Ri) , and include 
a new potential V 2 in a ring of internal radius i?i, and 
external radius R 2 that is used to tune the cloaking mech¬ 
anism (see Eig. [^. The basic idea is the following: we 
change V 2 in such a way that Sm vanishes, so that a ^ 0. 
Differently from a normal 2DEG, where it is necessary to 
tune Vi and V 2 to produce Jo = = 0, Eir [29], here it 
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Figure 2: Transport cross-section for a vacancy of radius 
Ri = 1 nm as a function E considering the 6 first partial 
waves (solid) and and a single partial-wave (dashed) . 


suffices to modify one of the phase-shifts, as Jq = are 
the main contributors to the cross-section; hence when 
they are zero, the total cross-section goes effectively to 
zero (it is reduced by several orders of magnitude). 



Figure 3: Sketch of the problem, showing the incoming 
electron, the impurity, and the cloak. 


We perform the same type of calculations described 
in the previous section but instead of defining the wave- 
functions in two regions, we now have three: Outside the 
potentials (region 3), for r > i? 2 , we have 


iA3J„+i(fc3r)e*(™+i)^ 

f Y^{ksr)e^^o 

For R 2 < r < Ri (region 2) we have 
( Jra{k2r)e^^0 

/ Y^{k2r)e^^^ 

whereas for r < Ri (region 1) we have 
, 1 . .. f Jm{kir)e^'^^ 


( 11 ) 


( 12 ) 


(13) 


where ki = \E — Vi\/hvF^ ^2 = \E — V 2 \/hvF and 
/c 3 = \E\/hvF^ The boundary conditions 7 / 1 ^(i? 2 ,^) = 


7 / 1 ^(i? 2 , 0 ) and '0^(i?i, 6>) = '0^(i?i, 6>) lead to four equa¬ 
tions. The solution of these boundary conditions deter¬ 
mine BjnlAjn and consequently 5m- 

To highlight the effects of cloaking, we choose Vi to 
maximize the initial cross-section, which is related to the 
resonant scattering through the disk of ratio Ri. Then, 
we tune the potential of the ring V 2 and its radius R 2 
in order to minimize the cross-section. Let us start with 
the situation described in the Figure where we have a 
disk with radius Ri = 10 nm and, for an incident beam 
with energy E = 0.02 eV (a value close to the Dirac 
point, which can be easily accessible for graphene on h- 
BN), resonant scattering occurs for some specific values 
of Vi. For Vi = 199 meV and Vi = —114 meV, the 
broad peaks in the transport cross-section are related to 
the resonance of the low-order phase-shift, (io = 5-i. On 
the other hand, for Vi = 282 meV and Vi =221 me, we 
can see a sharp peak in consequence of the resonance of 
the phase-shifts Si = S- 2 - Next, we tune the potential 
and the radius of the ring to investigate the maximal 
decrease in ctt in these particular situations. To identify 
the best cloaking scenarios, we analyze the behavior of 
the cloaking efficiency, which is the ratio between the 
transport cross-sections with, ctt, and without, ctto, tke 
ring, (Tt/o-to- 

Figure]^ a) shows a density plot of the cloaking effi- 
ciency as a function of the potential V 2 and radius R 2 
of the ring, for a fixed value of E and of the transport 
cross-section for the resonant scatterer uto = 132 nm 
(Fi = -114 meV). Panel gb), which is a lateral cut 
of Figure (a) for R 2 = 20 nm, reveals an impressive 
suppression of the transport cross-section of 3 orders of 
magnitude. 

In view of these results, we conclude that the cloak 
is very efficient in this setup. Here, we consider a sce¬ 
nario where both the resonant scattering and the cloak¬ 
ing are engineered with the use of tunable back and top 
gates [32]. In this particular case, it is necessary to work 
with disks and rings with radius in the order of tens of 
nm. Best cloaking efficiency, of the order of 10“^, can 
be obtained with thin layers with a distance of ^ 2 nm 
between internal and external radius. In spite of the fact 
that geometries with these sizes may be difficult to fab¬ 
ricate, fortunately there are some regions in the parame¬ 
ters’ space that provide a very effective cloaking for more 
realistic values of R 2 and V 2 - 

Moreover, for fixed i? 2 , there are large oscillations of 
the value of the transport cross-section as a function 
of the potential of the ring. These variations in the 
transport cross-section allow us to drastically change the 
transport regime: from cloaking to resonant scattering 
by applying small changes in the gate voltage V 2 - 

Figure(a) reports ctt/ctto as a function of the energy 
of the incident beam E and V 2 . We can see that cloak¬ 
ing is more efficient at low energies, close to the Dirac 
point, where the changes in transport cross-section are 
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Figure 4: (a) Cloaking efficiency as a function of V 2 and 
R 2 for E = 0.02 eV, Ri = 10 nm and Vi = —114 meV. 
(b) Cloaking efficiency for the same parameters of (a) 
and R 2 = 20 nm. 


Figure 5: (a) Cloaking efficiency as a function of V 2 and 
E for Ri = 10 nm, Vi = —114 meV and R 2 = 20 nm. 
(b) Cloaking efficiency as a function of E for the same 
parameters of (a) and V 2 = 286 meV. 


more pronounced. If now V 2 is set to minimize (JtIcftoi 
Fig. (b) shows the cloaking efficiency as a function of 
the energy and highlights the energy selective nature of 
cloaking: The original resonance was located ai E = 0.02 
eV, where the cloaking is more efficient. 

To obtain information on the direction of the scatter¬ 
ing process, we mapped the differential cross-section as 
a function of the scattering angle 0 and Vi (see Figure 

). The absence of back-scattering is very clear for any 
value of V 2 . Besides that, we can see that the cloaking 
is isotropic and there are changes in the scattering in¬ 
tensity according to the scattering angle in the resonant 
regime. We conclude that tuning V 2 allows the control 
of the scattering angle, demonstrating the feasibility of 
directional scattering. 

A. Cloaking vacancies 

We also investigate the cloaking mechanism for vacan¬ 
cies, a different source of resonant scattering. Again we 
keep the original vacancy with a radius of = 1 nm 


and include a new potential V 2 in a ring of radius Ri and 
R 2 {R 2 > Ri) , or equivalently, a disk of radius R 2 . We 
use the potential to tune cloaking, maintaining the basic 
idea: we change V 2 in such a way that it leads to cr ^ 0. 
The two boundary conditions = 'ip^{R2,0) 

and ^) = 0 lead to the ratio BmlAm, and conse¬ 

quently Sm- 

The density plot of the cloaking efficiency as a func¬ 
tion of V 2 and i? 2 , shown in Figure [^, reveals that the 
most efficient cloaking occurs for rings that are 1-5 nm 
thick. As discussed previously, thin cloaking layers based 
on back and top gates are difficult to engineer. Conse¬ 
quently, we focus the rest of our analysis on thicker rings, 
of the order of R 2 ^ 10-20 nm. For large values of i? 2 , 
there are strong oscillations of the scattering intensity, 
going from strong scattering to transparent states with 
small changes in V 2 ( see Figure [^). The distance be¬ 
tween the peaks in Figure is set by external radius R 2 . 
Remarkably, strongly scattering states closely coupled to 
invisible ones also exist in photonics systems, where the 
scattering cross-section of multi-layered plasmonic covers 
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Figure 6: Differential cross-section as a function of V 2 
and the scattering angle ^ for = 0.02 eV, Ri = 10 nm, 
Vi = —114 meV and R 2 = 20 nm 


may exhibit this so-called “comblike” behavior m- Here 
a similar behavior in the cross-section occurs as a func¬ 
tion of the external potential, and not as a function of 
the incident wavelength as in the photonics case. This 
result would allow one to dramatically modify electron 
scattering in graphene by varying an external parameter. 

An important difference between the disk and the void 
cases is the selectivity in energy. In the first scenario, 
there is a reduction in the value of cft of 3 orders of 
magnitude at the specific energy where the resonant state 
is located, which in turn can be tuned by the potential 
Vi. For the vacancy, the resonance is always located at 
the Dirac point and cloaking only occurs in the vicinity 
of = 0 (see Figure a). For energies very close to 
the Dirac point, the cloaking is very efficient, as shown 
in Figure [^b. However, for the energies discussed in the 
previous section, the transport cross-section is reduced 
to only 10% of its original value. 

Considering just the first partial-wave. It is possible 
to obtain a simplified expression for maximum cloaking 
(^0 = 0) at the Dirac point (F^ = 0). The condition reads: 

Ji{x)Yo{ax) — Jo{ax)Yi{x) = 0, (14) 

where x = and a = ^. This equation is satisfied 

for a sequence of increasing values of x for a given a < 1. 
In the limit x ^ 1 the above equation acquires the form 


cos[x(l — (a)] « 0 , (15) 


whose solution is 

(2n + l)7r 
2(1-a) 


(16) 


with n = 0,1,2 ,..., and it is consistent with Figures 
andlHl 




V2 (eV) 


Figure 7: (a) Cloaking efficiency as a function of V 2 and 
R 2 for E = 2 meV and a vacancy of radius Ri = 1 nm; 
(b) Cloaking efficiency for the same parameters of (a), 
F^ = 1 meV and R 2 = 20 nm. 


B. Condition for cloaking by an arbitrary potential 

In this subsection an integral equation for the phase- 
shifts is derived and a general condition for the vanishing 
of the phase-shifts is obtained. For a general potential 
V (r) we seek solutions of the Dirac equation in the form 

i’m{r,e) = ) • ( 17 ) 

Acting with the Hamiltonian Rq -j- V (where Rq is the 
Dirac Hamiltonian) in the above general form of '0m(^5 ^) 
we end up with two coupled first-order differential equa¬ 
tions: 


G'+'^^^G = eF, 
r 

(18) 

F' -—F = -eG, 
r 

(19) 

where e = {E — V) / {xph). We benefit from symmetrizing 

the two last equations by introducing two 
G = g/\/r and E = f /y/r. These lead to 

new functions. 

9' + -9 = ef, 
r 

(20) 

H 3 n 

f - - f = -^ 9 - 

r 

(21) 
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Figure 8: (a) Cloaking efficiency as a function of V 2 and 
E for a vacancy of radius i?i = 10 nm and R 2 = 20 nm. 
(b) Cloaking efficiency as a function of E for the same 
parameters of (a) and V 2 = 87 meV. 


where j = m + 1/2. For zero potential, we write 

g'o + ^go = eo/o , (22) 

/o - bo = -eogo , (23) 

where eo = E/{vfK). We can now combine the two sets 
of equations, with and without the potential, by an ap¬ 
propriate multiplication of the equations with the inter¬ 
action by the spinor components with zero potential and 
vice-versa. This procedure leads to 

^(/o5) = /o/e-55oeo , )|;(5o/) = -goge+ffo^o , (24) 
which we can subtract from each other to obtain 

■^ifog - fgo) =-:^iggo +fof) ■ (25) 

Integrating the preceding equation from zero to infinity 
and recalling that /o and go vanish at r = 0 we obtain 

lim [fog-fgoT° = -b f drV{r){ggo+fof) (26) 

a^o^oo Vpa Jq 


We also note that for r ^ oo, the scattering states are 
simply related to the asymptotic forms of the Bessel func¬ 
tions as 


(27) 

(28) 


f ^ \l cos{kr - Xm + dm) , 


^ ^ \l sm{kr - Xm + dm) ■ 


Using the latter asymptotic forms, Eq. 
sin dm = -r 

/o 


reads 


kTT 


1 f°° 

sin dm = -r / drV (r) {ggo + fof) ■ 

Jo 


(29) 


Now /o = x/rJm{kr) and go = X^/rJmFi{kr), and we 
obtain 


I- ^m, — 


nk 1 
2 vph 


nOO 

/ dr^V{r)[fJm{kr)+gXJm+i{kr)]. 

Jo 

(30) 


We note that there is an alternative route for deriving 
the last relation, but it is a much more complex pro¬ 
cedure that uses the Lippamnn-Schwinger equation and 
the Green’s function of the Dirac Hamiltonian. For a 
potential of finite range r = i?, we have 


sin 6 rn 


irk 1 
2 vpk 


nR 

/ drV^V{r)[fJm{kr)+gXJm+i{kr)]. 
Jo 

(31) 

From this equation we extract the general condition that 
is sin(5m = 0 {cosdm = 1): 

p 

0 = drVrV{r)[fJm{kr)+gXJm+i{kr)]= Km, (32) 
Jo 

which is the central result of this subsection, since it cor¬ 
responds to the general invisibility condition for an arbi¬ 
trary potential V (r). Let us consider a simple application 
of the above result taking into account a potential that 
is constant up to r = R. Then Eq. (32) leads to 




drr[Jm{pr)Jm{kr) + JmFi{pr)JmFi{kr)] = 0 , (33) 


where p = {E — V)/{vfK) (we assumed E > 0 and E > 
V). Performing the integral, we obtain the condition 

0 — pJjYi—\ {pRJj JjYi{kRJ) kJjYi—\ {kp^Jmip '^) 1 ( 84 ) 

which is the same we obtain from the phase-shift analysis. 


In the above form, the result (32) is both a condition 


for the vanishing of the phase-shift caused by a potential 
and a cloaking condition if the potential can be sepa¬ 
rated into two (or more) different contributions that can 
be tuned separately, as discussed in the previous exam¬ 


ples. Finally it should be appreciated that condition (32) 


depends only on the value of the wave function in the re¬ 
gion where the potential is finite. In the next subsection 
we discuss an intriguing problem, the ability of a spatially 
short-range potential to cloak a resonant scatterer. 










C. Cloaking a void with a (5— function potential 

Let us consider the cloaking of a void by a (5—function 
potential of the form V{r) = VoS{r — Rq); we also de¬ 
fine vq = VoKvfK). The fact that the Dirac equation is 
a first-order differential equation and the (i—potential is 
rather singular implies that the boundary condition can¬ 
not be the continuity of the spinor at the location of the 
potential. The radial part of the Dirac equation can be 
written in matrix form as 

d<S 

= —iecFyS -f icFyvS — -cfzS = G{r)S , (35) 

where = (f^g) and v = V/{vfK)- Formally, the latter 
equation can be integrated as 

S = rre^^^o^^'^(^'^S(ro), (36) 

where % is the radial position-ordering operator, defined 
as 


7;[G(ri)G(r2)] = 0(ri - r2)G{ri)G{r2) 

+ 0{r2 - ri)G{r2)G{ri). (37) 

Performing the integral, which is now a well defined quan¬ 
tity, the spinor right after the position of the function 
reads 




S{Ro^r]) = TrC^^o-ri 


S{Ro-v)=e^^y^^S{Ro-v), 

(38) 

and 7^ = 0+. Eq. (38) embodies the boundary condition 
obeyed by the spinor at a (i—function potential. Using 
Lagrange-Silvester formula, we can compute the value 
of the exponential as follows: for a function f{A4) [one 
should not confuse the spinor component / with f{M)] 
of a 2 X 2 matrix Ai we have that 


f{M) = /(mi) 


1712 - M 

1712 — 


+ f{'rn2) 


mi — M 

mi — m2 ’ 


(39) 


where mi /2 are the two eigenvalues of the matrix A4. 
Applying this result to our case we obtain 

f{iayVo) = cosvq + iay sinuo . (40) 


analysis. We should stress that for this problem the “free” 
solutions are those referring to the scattering states in the 
presence of the void. Therefore, Jm{kr) and Jm-\-i{kr) in 
Eq. (30) should be replaced by the spinor components 


of the scattering states in the presence of the resonant 
scatterer alone, that is, by: 

Jm{kr) cos5^[Jm{kr)+ta.nSl^Ymikr )], (42) 


where tan 5^ is defined by Eq. (10). In the presence of 
the (5—function potential, the spinor to the right of the 
potential has to have the form 


'>Pm{r,0) = cos Sr, 


Jm{kr)i 


A9m 


— tan(5„ 


(43) 


iJm+i(A:r)e*^('"+i) 
Ym{kr)e^^^ 

At low energies, we find the following conditions that 
make the first two first phase-shifts vanish: 

r cosuo — + 2 \n{R/Ro)] sinuo = 0 , / = 0 

I — cosuo + |i?o(3i?^ — 2i?Q)/c^ sinuo = 0 , 1 = 1 

(44) 

Note that for vq ^ 7r/2 the cloaking condition is purely 
geometric: for (ig = 0, we have the condition Rq = ^/eR, 
whereas for = 0 we obtain the relation Rq = y^372i?. 
The two numbers do not differ much from each other. 
Therefore, we can almost cancel the two phase shifts at 
the same time. 

In Eigs. an d pT| we depict density plots of the func¬ 
tions \Kq\ and \Ki\. The regions where \Km\ = 0 corre¬ 
spond to the vanishing of sin Sm and therefore to a cloak¬ 
ing effect on the void by the (i—function potential. In 
Fig. |9]a density plot is depicted as function of Rq and 
k. The vertical dark blue regions with \Km = 0| corre¬ 
sponds to the Rq = y/eR and Rq = y/3/2R for \Ko\ and 
1^11, respectively. Interesting enough we see that these 
geometric conditions for cloaking hold as long as kR < 1. 
In Eig. a density plot is depicted as function of Rq 
and Vq . As in the previous figure we see dark blue regions 
where \Km\ = 0 around vq « 7r/2. 


IV. TIGHT-BINDING CALCULATIONS 


Thus, the boundary condition obeyed by the spinor reads 


S{Ro + 7/) = (cos7;o + icTy smvo)S{Ro - rj ). (41) 

We note that in this case the spinor is not continuous, as 
we had anticipated, and therefore we need a prescription 


for integrating the function in Eq. (30). To that end. 


we have to assume that both spinors, that to the left 
and that to the right of the potential, contribute half to 
the integration. Eurthermore, we should not forget that 
the spinors at the right and at the left of the function 


are connected by Eq. (41). If this is done, then the 


previous method yields the same result as the phase-shift 


To confirm our predictions and to test the robustness 
of the cloaking against disorder, we also model our sys¬ 
tem with a tight-binding Hamiltonian, and use the ker¬ 
nel polynomial method (KPM) [44] to calculate the local 
density of states of a graphene sheet [45|. This real-space 
approach allows us to introduce gated regions with arbi¬ 
trary geometries. Thus, we consider the possibility of 
having a gate potential without perfect radial symme¬ 
try by introducing roughness to the ring potential. The 
Hamiltonian reads 

n = + U ^ cla + V 2 J 2 4g (45) 

(i,j) iEiRi iEiR 2 
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Ro (A) Rg (A) 


Figure 9: Density plot of \Krn\, with R = 2.8 A, and Figure 10: Density plot of \Krn\, with R = 2.8 A, and 

'^0 = 7r/2, as function of k and Rq. Top panel: lA'ol; k = 0.034 1/A (corresponding to V/ = 50 V in a 

bottom panel: \Ki\. graphene FET on Si02), as function of vq and Rq. Top 

panel: | Kq \; bottom panel: \Ki\. 


where q is the annihilation operator of electrons on site 
i where t ^ 2.8 eV is the hopping energy between nearest 
neighbors (NN) sites in a honeycomb lattice. Vi is the 
radial potential and V 2 is the potential of the ring. 

For a tight-binding hamiltonian with N=2 x 600 x 600 
sites and periodic boundary conditions, we follow the 
same scheme that was discussed in the previous sections: 
we consider a disk of radius Ri = 8a where a is the lat¬ 
tice constant of graphene and tune Vi to obtain a quasi¬ 
bound state inside the disk, as illustrated in Figure pT] for 
Vi = 0.45t, which is equivalent to the resonant scatter¬ 
ing regime. We then vary V 2 inside the ring with radius 
Ri = 8a and Ri = 16a to maximize the cloaking. In Fig¬ 
ure E] b, the quasi-bound state disappears, highlighting 
the cloaking efficiency (even in the presence of disorder, 
as we shall see in in the following). Finally, we introduce 
roughness in the external side of the ring and compare 
the local density of states with (Figure pT]c) and without 
disorder (Figure pT]b). The cloaking is highly efficient in 
both situations, demonstrating the robustness of our re¬ 
sults against disorder, which is intrinsic to the fabrication 


process. The same procedure can be repeated in the case 
of the vacancy and the results were found to be similar. 

V. CONCLUSIONS 

In conclusion, we have envisaged novel, alternative 
electronic cloaking mechanisms in graphene. These 
mechanisms are based on drastically reducing the scat¬ 
tering cross-section of resonant scatterers in graphene, 
which can be either an engineered resonant state pro¬ 
duced by a top gate or a vacancy. Using a partial-wave 
expansion of the electronic wave-functions in the contin¬ 
uum approximation, described by the Dirac equation, we 
demonstrate that electronic cloaking is very efficient; the 
reduction in the scattering cross-section can be as im¬ 
pressive as 3 orders of magnitude. In the particular case 
of the resonant state, we show that cloaking is strongly 
facilitated by the fact that the first phase shifts, which 
give the dominant contribution to electron scattering, are 
identical for graphene. Hence one just need to tune one of 
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Figure 11: Logarithm of the local density of states for 
E = 0.00 eV. Panel (a): Vi = OAbt V 2 = 0. Panels (b) 
and (c) Vi = OAbt V 2 = 0.23t. The sites with Vi = Vi 
are black while the sites with Vi = V 2 are green. 


the scattering potentials to reduce significantly the scat¬ 
tering cross-section, this in contrast to other electronic 
cloaking schemes in graphene proposed so far [29]. We 
also show that sharp transitions between invisible and 
strongly scattering states occur for some values of the po¬ 
tential. We suggest that this result, which can be viewed 
as an electronic analogue of optical comblike states in 
coated scatterers [31], could be explored in applications 
aiming at tuning the electron flow in graphene. Using 
a symmetrized version of the massless Dirac equation, 
we derive a general condition for the cloaking of a scat- 
terer by a potential with radial symmetry. By means 
of tight-binding calculations, we show that our findings 


are robust against the presence of disorder. Altogether, 
our results reveal that electron cloaks are not only more 
easily achieved in graphene, but also that they could be 
applied to selectively tune the direction of electron flow 
in graphene by applying external gate potentials. As a 
result, we hope that could serve as the basis for designing 
novel, disruptive electronic devices. 
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